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ABSTRACT 


Electron impact ionization cross section, energy distribution and the 
average energy of the ejected electrons for the first five principal quantum 

numbers of atomic hydrogen, and all the azimuthal quantum numbers are cal- 

% 

culated within the Born approximation. Calculations are extended to a range 
of 100 to 900 threshold units of energy. Regions of agreement and disagree- 
ment between the high energy quantal and classical calculations are analyzed. 
Validity of the Born approximation for excitation and ionization for a single 
particle model is given. An energy spectrum of the secondary electrons due 
to the passage of an energetic electron in a hydrogen atom gas is found. 
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X. INTRODUCTION 

Ionization of the ground and excited states of atomic hydrogen by electron 
collison has a number of applications in plasma physics and astrophysical prob- 
lems. While there has been an extensive study for the ionization of the ground 
state, little has been done for the excited states. Among studies for the ground 
state ionization we can mentionthe early calculations by Bethe (1930), and 
Massey and Mohr (1933) within the Born approximation; measurement of the 
cross section by Fite and Brackman (1958), and Boksenberg (1960); measurement 
of the cross section at the threshold by McGowan et al. (1968); study of the wave 
function for e-H ionization by Peterkop (1962); distorted-wave calculation by 
Veldre and Vinkalns (1963), and Burke and Taylor (1965); and study of the threshold 
of ionization by Wannier (1953), Geltman (1956), Budge and Seaton (1965), and 
Temkin (1966). 
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Here we are concerned with ionization of the excited states. In a previous 
paper (Omidvar 1965) a calculation on impact ionization of the excited atomic 
hydrogen was reported in which the cross section and other related values were 
given for wave functions in parabolic coordinates. In the present paper these 
values are given for wave functions in spherical coordinates. In problems re- 
lated to radiation, these are more suitable coordinates. Similarly, in more 
complicated calculations than the Born approximation, spherical coordinates 
are more often used. For these reasons the present calculation will have a 
wider range of applicability. 


In the present paper, calculation is extended to much higher incident energies 
than used in the previous paper. For principal quantum numbers n ~ 1, 2, 3, 
calculation has been extended to 900 threshold units; for n « 4, to 225 threshold 
units; and for n = 5, to 100 units. The numerical results presented here are 
more accurate than the previous calculation. Finally, in this paper, energies 
are expressed in threshold units, and scaling law has been used for other 
quantities. Thus it has been possible to study more precisely the variation in 
cross section and other physical quantities in terms of n and the azimuthal 
quantum number / . 

On the other hand, values for different magnetic quantum numbers m are 
not reported here. In this paper the axis of spatial quantization is taken along 
the momentum of the ejected electron, while in the previous paper this axis was 
along the momentum transfer of the incident electron. Values for a given m 
in the two cases naturally do not correspond. Unless one is interested in the 
ionization cross section and related values for a particular momentum transfer, 
or for a particular direction of the ejected electron, the m dependent values 
serve no useful purpose. Values related to a given n and / are listed here 
after an integration is carried out with respect to both the direction of the ejected 
electron and the direction of the momentum transfer. 

An analysis is also made for the validity of the Born approximation and 
its relation to the Bethe approximation. Also a comparison is made between 


high-energy quantal and classical calculations, and it is shown that as n increases 
in accordance with the correspondence principle, the quantal results approach 
those given by the classical calculation. 

2. DERIVATION OF THE IONIZATION AMPLITUF 

The analysis must start with the evaluation of the Born amplitude for 
ionization for any given state. For this it is necessary to specify the states of 
the atom before and after ionization. Atomic units are used in which if and the 
electronic mass and charge are set equal to unity, and energy is expressed 
throughout in rydbergs. If p and e specify the momentum and the energy of a 
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particle, and k, its propagation vector, it then follows that p ~ k and e - k . 

In this way energy takes the dimension of inverse of the length squared, but 
numerically its value is given in rydbergs. 

States of the atom are specified by the Coulomb functions. The initial 
state is given by \nlm> , where n, / , and m are the principal, the azimuthal, 
and the absolute value of the magnetic quantum numbers. This state is related 

to the states in the parabolic coordinates, \nn l m> , by the unitary transformation 

n - m - 1 

I nlm > » 2 I wi, m > < nn m | nlm> , m 

n x « 0 
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where <nnm | nlm > are elements of the transformation matrix. In this 
transformation n and m remain fixed, «, ranges from 0 to n-m -1, and / 
ranges from m to «-l. The recursion relation and the specific values of 
< nn { m | nlm > for n = 0 to 5 and all possible values of m are given elsewhere 
(Omidvar, 1967). Barut and Kelinert (1967), and Hughes (1967) have shown 
independently that these elements are related to the Clebsch-Gordan coefficients 
through 

m 1/2 

< nnm | nlm > - (-) (2/ + 1) 


/l/2 («-l) 1/2 (w-1) I 

\l /2 (m-n , +« 2 ) 1 12 (w+« , -« 2 ) *w 



m ’ n x +n 2 m+n \' n 2 


2 


2 


/ /, / \ 

where I I are Wigner's 3/, and C(i ji; mm„m) 

/ 12 1 2 

Clebsch-Gordan, coefficients (Edmonds 1957)*. 


mj ( 2 ) 

are the 


In (2) the integer is given by n 2 s h -1-m- n ( . The normalized hydrogenic 
state | nn^m> corresponding to a nuclear charge Z is given by Bethe and 
Salpeter (1967) 


♦Equation (2) is given by Barut et al. If we write A[ n for < nn^m \nlm>, then it 

can be shown directly that A f = (-) n 2 A ^ » O -)” 2 * ^4 . , where 

in j in j in | 

ln x ’ /«, 


i ’ i 

, and are the matrices of Barut et al. , Hughes, and this 


author, respectively. These phase differences have no effect on the calculation. 
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In this equation 1 ?, v? » are the parabolic coordinates of the atomic electron 
related to the spherical coordinates by £=r (I + cos 0 ) , rpr ( I - cos 0) and ^ ; 

and a„ is the Bohr radius. 


The final atomic state is the state of an ejected electron in the field of a 
nucleus of charge Z, and is given by a Coulomb function with outgoing waves. 

If k represents the propagation vector of the ejected electron and at the same 
time the direction of the z -axis of the coordinate system, this function is given 
by (Landau and Lifshitz 1958) 


I k > = ♦(, + (r).,V c c ' k ' 1 F W, I, ikv) , 


^ c " 2* [/ - exp r-2*fij 


,m k- 


where F(a, b, z ) is the confluent hypergeometric function which is regular 

+ 

at the origin. ^ (r) is normalized such that 


\|£<r) tp* (r) dr. S (k-k) - 2k' S(e-e’) «(£-£') , 
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with e and e' the energy of the ejected electrons in rydbergs. The asymptotic 

Hh 

form of $ (r) is given by 


* + (r) 


(2n) 


•3/2 ikr 


(6) 


In some application it is more useful to use an integral representation for \p (r) 
which originally was given by Sommer feld (1931). This is 


+ e 

f. (f)-N c - 

k p 


s ikr f°° , u .ip ___ 

I e u J a (2y/iukr | 

u-mi 


) dU , 


(7) 


However, in this paper we use the form given by (5). 


With the atomic states so defined and with k t and k 1 the propagation 
vectors of the incident and the scattered electrons, the ionization cross section 
per unit energy range of the ejected electrons is given by 


dQ m !_[ 

de 4ira 2 Q k 2 I 


kd k 


A + k 2 

\ -k 7 


I T (/, f)\ qdq, 


( 8 ) 


where the transition amplitude TUJ) is defined by 


raj) 


4n 


<r 1 1'*' ■ r i /> . 


( 9 ) 


with q - kj - k 2 5 and the total cross section defined by 


fmctx 

I ^ 

I de 

J a 


de 


(10) 


The values of k 2 » e , and € max are related to the incident energy k 2 by conser- 
vation of energy through 


A?* -Z 2 / n 2 a J - k* + e = e 


max 


(H) 


with n being the principal quantum number of the initial state. To have a clear 

account of the dimensions involved notice that with reference to (3) and (4), 
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T(i,f) has the dimension of HI . Then dQ/de and Q have the dimensions 
of the fourth and the second powers of length respectively. 

The rest of formulation consists of evaluation of T(i, f) given by (9). 

We take! / > and |/> to be the states given by (3) and (4). For simplification 
in evaluation of T ( i, f) , we take the z -axis of the coordinate System not 

along k^, as is customary, but along k. We represent the spherical coordinates 
of q in this system, by q, 0 , ^ and the parabolic coordinates of r by q, <p, . 
The ^-dependent integral in (9) and its value ate then given by 
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2 7T 


( iq • r + i m *p) , 
e dip 


-2n exp ft i q (cos 6 1 ) ($ -q) + /m^ ( 4 Vi it) J m ( q sin 0 J \f%r\) , 


( 12 ) 
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whore (*) is the Bessel’s function with argument * . Using this 
expression we c an write 


7 (mi. m, k)* 
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i q * 


2s - - ,m( ^ J) N n „ m N c — 

"I 11 !™ C d/A 
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p&Aot , s - -‘A l(q cosO x -k) , 


(14) 


To evaluate (14) we use the identity 


r M/ ^ H m n\ z „ m d"-™ 

V* *4^' * 




Then by direct differentiation 


m AAA 

I (or*) » (-) . 

Aij + m 


n i f 


(15) 
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Vj = 0\ 1/ (r | -f AW )/ 

With this form for L m (ot%) , the (--integral in (14) can be evaluated and 
its value is given by 

/ \ v-. ■ "•/{ v 
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m + J 
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r P* p+s. 


msasm 




When this is substituted for the £ -integral in (14) we obtain 


Bm( f (», +>»)< II f w -> ,) 

~rr~Zj\ p J + "*» <u*> \ (2 />)"'+ 1 ) 

1 p r ° 

/ - / dn 1*^ i\ n i 1 (an) F (•$> 1> dkn ) 1 > 

JO 1 n 2* m J 

\ - P‘S + q sin 9 ^4P . 


Evaluation of the integral in (17) is accomplished by expanding L (an) 

n + m 

2 

according to (16). The integral can then be evaluated by standard methods 


(Landau and Lifshitz 1958), The result is 


m + v 


2 „”o V V ' ! 


(« 2 + m)! • /«A (« 2 + m)! v /m + Vj 

1 ~ nj (•',) (<\ + »l)! ” * \ M, 


«. = ° 


y (-V-D ! OP)! y (X + ikf-Pj 

(-<0-1 -m-« 2 + « 2 )! (/0- m 2 )! A x /0 + 1 + m + « 2 - m 2 - 



By some manipulation and using the properties of the binomial coefficients, this 
can be written 



Substitution of / in (17) and evaluation of U and bU/bp is straight forward 


algebra. By inserting the value of bU/bp in (13) we find that 
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with 


-M - v -n - v n -M t + M 2 + ^ ^ 

F'/i 4 ‘ 2? D D* 1 2 1 
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2{ct + /*) (m+M 4 +V, ■ W) 2«(i +m 2 +^ + //3) M 2 ' M 4 -M,+M 4 +^+^ 

• * — + — + — — 

A B D D* 


where D* is the complex conjugate of D . In (22) the lower limits of the 
variable integers are zero, and in (23) A, B, and D are defined as 

A = (Q' + /jt) 2 + q 2 , B*at 2 +( q-k) 2 ,/)= a+£»(q-k) , 


In this way, in the expression for T (n n m, k) all symbols are defined except 
C( ji .n ). This coefficient arises from the following consideration. 

nr 


Suppose 


/IP 

U (P) = 

-PP + C 


where A, B, and C are constants. Then the l derivative of y = y (u) with 
respect to P is given by 






* ( 23 ) 


(24) 
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(26) 


fly l 

-p-Ws) 1 S c(\n 

X* 1 


AC A X 


(¥*) 


^ X («), g- 


•i/P+C 


I 


X //• 

where y A (w) is the X' derivative of y with respect to u, and C(\l) are numbers 
depending on X and /, Differentiating equation (26) on one hand with respect to P, 
and, on the other hand, substituting / + 1 for l in this equation and equating the 
resulting expressions, we find the recursion relation for C(X, /) is 


C(X,/)» (M+X) C (X, /- 1 ) + C(X-1,M), X</. (27) 

We also note from (26) that C ( 1, 1 ) = 1 ; and we put C(X, / ) = 0 for X>/. These 
conditions and (27) are sufficient to determine all C(X, /). 

It is convenient to define C(0„l ) • 5 (0,1). Then (26) can be written in a , 
more general form 

d l y , 1 /AC \ X X 

~P‘ (Bg) 1 £ C(X,/) g\ y (£/), (26a) 


with g defined as before. The use of C(X, /) becomes necessary when in (17) 
differentiation with respect to P is carried out to evaluate £/. A • numerical table 
for C(X,/) is given elsewhere (Omidvar 1965). 



The expression given here for T(nn m, k), is also given elsewhere 
(Omidvar and Sullivan, 1967)*. 

Since the transition amplitude is linear with respect to the initial state 
wave function, from (1) the transition amplitude in spherical coordinates, 
T(nlm, k) , is given by 

n ~m-l 

T (nlm, k) - £ <nn l m | nlm> T (wyn, k) . (28) 

= o 

This completes the evaluation of the transition matrices. 


3. VALIDITY OF THE BORN APPROXIMATION AND BETHE APPROXIMATION 
In collision of a charged particle with an atomic system, the use of the first 
Born approximation amounts to two separate approximations. One is the assump- 
tion that the atomic system has only one mode of excitation or ionization. The 
second assumption is that the incident particle moves in a force-free field and 
is not affected by the atomic potential. Physically there are infinite modes of 
atomic excitations or ionization, virtual or real, depending on the energy of the 
incident particle. The validity of the first assumption then depends on how 


*The phase exp [im (<p + Vxn)] does not multiply T (nn m, k) in this 
reference. This phase becomes important when we consider the effect of the 
electron exchange. 
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strongly the excited mode in question is coupled to the neglected modes of 
excitation. Mathematically, its validity depends on the convergence, and the 
rapidity of the convergence, of the higher order of the Born approximation 
(Moiseiwitsch and Smith, 1968). 

The second approximation is analoguous to the scattering of a particle by 
a single center of force, and similar to the case of the single center of force , 
a criterion for its validity can be derived. This will be done here. 


By designating the position vectors of the bound and the incident electrons 
by r j and r 2 , the exact transition matrix will be given by 


r E l 1,2). (r 1 )e'' S ' r * 


1 Z 


r, „ r n 

I 2 2 


* (t,,r 2 )> 


where 1 and 2 refer to the initial and final states, (r J# r 2 ) is the total wave 
function of the system, and <P 2 (f ) is the final state atomic wave function. 
The two approximations mentioned previously consist of the replacement of 


* (r i’ r 2 ) 


by ( r j ) exp(tk ,r 2 ) . Thus when integration 


is carried out with respect to iv , (9^ ensues. However, without carrying out 
the integration, (29), with this replacement, can be written 


r i «i • r 

r (1, 2) = J V (t) e dr 


v <r 2 ) - <M r ,> 


r - r 

1 1 2 


0! (Tj) >» 
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where the Z / r 1 is dropped on account of the orthogonality of 0, and 4> 2 . (30) 
is the scattering amplitude of a potential V (r) in the Born approximation 
(Schiff 1955). Corresponding to this amplitude the scattered-electron wave 
function, u (r) , is given by 

«(r). e‘ k ' ' r + W (r) (32) 

where 



/k r/ + //r | r-r ' | 
g 1 — 1 dr > 

Ir-rM 


(33) 


The condition for the validity of the first Bom approximation is that | w( r)|< 1 
for all values of r. 

In the problem under consideration, besides the short range forces, an important 
contribution arises from the long range force due to the 1 /r 2 potential. The 
following case considers a bound-to-bound transition, and derives a general 
expression for V {t ) when r is small or large, Then, by making use of (33), 
the validity of the Born approximation is established in terms of the energy of 
the incident particle. The results obtained can also be extended to the case of 
ionization. 


. In spherical coordinates, the discrete hydrogenic-wave function is given 


*(»/»!, r) - r 1 P (nl. r) Y , (A), 


where Y, (P) is the spherical harmonics, and r H P(nl,r ) the radial 

1 m 

function given by (Landau and Lifshitz 1958) 


P (nl, r) * b e-tip pl+l F(-n+l+ 1, 21 + 2,9), 


i L r a (n + 1) i 1 

c 21 + 1)1 I2n(n-I-J)\] 


2Z 

a m •- — , p = ctr 

na 0 


For evaluation of (31) we make use of the standard expansion 


r i * r 2 


■ E E — ^ r. (A) y* (A) , 

>.--0 + l X " 


with r< the smaller, and /*> the larger, of jr^ and r 2 . Substitution of (34) 
and (36) into (31), and an integration with respect to P t leads to 

\ 4 * < 2 /,+/ ) 1 ^ 

s ? [(2/ + ;)(2X + 7)J • C (,/ H m W) 


C (/'X/, m'pm) 


X cm 000 ) y^ (P 2 ) w („/«'/' r,) , 


where 9= m- m' , and \ takes on values for which the C coefficients do not 


vanish. Also 


17 


'I 

0 nln'l ' ; r 2 ) m rf^ + D J P (nl, r 2 ) P (n'l f , r 2 ) r 2 x dr 2 

0 

OO 

+ r if P Oil r t ) e (n'l'.rj r’^ + 1) , (38) 

r , 

First we consider the case r < 1, Since in this case P (nl, r ) varies 

as r 1 + * , and always X < / + 1* , (38) leads to 

OO 

U x (nln'l', r<l) ~ r X f P (nl, r') P in'l', r') r> + 7) dr ' , (39) 

* o 

The definite integral on the right-hand side can be evaluated for specific cases. 
Since the smallest value of X is \l-l' \ , we conclude that for a bound -to -bound 
transition the potential is finite at the origin with the proportionality relation 

V (nln'l 1 , r) « r I I , r+o . (40) 

In contrast to the case of elastic scattering where the potential is a screened 
Coulomb potential, the inelastic collision potential is finite at the origin. 

For the opposite case of r> 1 we obtain 

OO 

1 ^^(nln'l', 1) ^ P (nl, r 1 ) P (n'V ’ r 1 ) dr' , (41) 

Since the definite integral is a constant, and the smallest value of X is | l-V | 
through (37) we find 

V (nln'l', r) <x r m \ l m l f \ m l t r + oo, (42) 
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"y* to-**** 
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An exception is the / » V transition in which case the integral in (41) vanishes 
on account of orthogonality of the radical function, and the potential vanishes 
exponentially. The potential in this case is then similar to a square well 
potential. 


From (42) we conclude that when / and l' differ by unity, the incident 
particle, in addition to short range forces, is also moving in a J/r 2 potential. 

An expression is drived below for the validity of the Born approximation for a 
dipole potential o* /r 2 , with V being the dipole moment. But first cr is evaluated 
for transition between the ground and any excited state with M (primes are 
dropped for the final states for convenience). Evaluation of cr follows from sub- 
stitution of P(1Q, r) and P(nl, r) from equation (35) into the integral in (41), 
and a straight forward integration. By taking Z* 1 we find that 


VilOOn Jm, | r | > 1) 


y* 

3 1m 


Y* (r) pr 


16 n 

cr " nr,.* Z. s /2 


s[f\ {n 2 - 1) \n + 1 


Values of cr in units of e 2 a Q for n= 2 to 10 are given in Table 1. From 
(43) it is seen that as n increases, cr diminishes as n" 3 / 2 . 


Considering the case of ionization the dipole moment depends on the velocity 

of the ejected electron. To treat this case, first a partial wave expansion of 

. + : 

& (r) , (4), with the same normalization and phase, is made. With the z-axis 

along the momentum of the ejected electron, this is given by (Schiff 1955) 



r 




* ( kl, , r) - /"> £ P m, r) Y, 0 (P) , 
1-0 


where 


p (kl, r) * Cfci e~^P + J F (-1(3 + 1 + 1 , 2 1 + 2, p), 


N ( c J4¥(2T+ 1) rrj , Z 

** ▼ II /irt \\ if. /■* 


a (21 + 1 ) ! 


JJ (IP -X), a - -2/X', 0 


,P - 


X-/ 


with jV c given by (4), It should be noted that by letting i(3-*n , a will become 
the same as that given by (35). 


To evaluate the dipole moment for the case of ionization, P ( nl r) in (41) 
is replaced by P(10, r) given by (35), and PQi'P, r) is replaced by P(kl, r) 
given by (45). Then for X = 1 the integration is carried out by the standard 
methods. In this way by combining (37) and (41) it can be seen that* 


V (lOOklO, \ r\> 1) ». cr cos (r . k) fr % , 


8 s/2 I 2n/k 


(1 + iM exp [ -2ip/k) 


2tt \l -exp [-2n/k] f (1 + k 2 ) 3 
For k< l, v - 8€* j\pnk. , and for k<l,v*i4 s/jj/vk * . 


", ipmtatf* k 


*When an electron is scattered by a potential which has imaginary parts, 
conservation of the current for the ingoing and outgoing wave functions is not 
satisfied. The lack of this conservation law also holds for a real potential in 
the Born approximation. 
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To summarize, in an inelastic collision of a charged particle with an 
atomic system, the effective potential in which the particle moves is always 
finite at the origin. Except for /* l' transitions, the potential has long ranges, 
the most important being Iff 2 for /'■ / 1 1 transitions* 

Since the potential is finite at the origin for small distances from the 
nucleus, it can be a*# ^oximated by a square well potential. The criterion for 
the validity of the Born approximation for this potential is that k l a> 1, 
with the momentum of the incident particle, and a the potential range 
(Schiff 1955). It will be shown below that the criterion for the validity of the 
Born approximation for a potential 



V{r) 



r>a 

r 


(47) 


is similarly k x a> 1 , 

To find the criterion for the validity of the Born approximation for the 
scattering of a particle in a — potential, we replace V (r) in (33) by its dipole 
form given by (41). By introducing x * r-r% this leads to 


# 




(48) 
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For the purpose of angular integration in (48) we take the z» axis along fc . Then 
making use of expansion (36) and the expansion 


OO \ 



and Lifshitz 1958). This will give rise to resonances not predicted in the Born 
approximation. The condition k x a> 1 also insures that this resonance will not 
happen. 

For potentials falling faster than V r 2 , the criterion k { a> 1 is naturally 
sufficient. Then by talcing the effective potential as the sum of the long and short 
range forces, with the neglect of the other modes of excitations, the Born approxi- 
mation is valid when k^aP 1 . 

The range a for different transitions can be taken to be the distance 
beyond which the product P(nl, r) P (V/ 7 , r) falls off exponentially. With reference 
to (35) we may take 

a • 2 (a + ex')' 1 . (53) 

For excitation this condition reads a = ( n + n 1 ) * nn'(ajZ) . Then for 
excitation of the ground state it is found that aszajz , while for excitation of the 
highly excited states, when n'=n, it is found that a = n(a Q /Z ) . In the 

latter case the ranges expands. The condition k P 1/a then implies that for a 
fixed incident energy, the Born approximation is more valid for excitation of the 
excited states. In the excitation of the excited states it should also be noted that 
the threshold energy is smaller and the peak in the cross section appears at 
smaller incident energy. If we measure the wave number of the incident electron 
in units of the wave number of the bound electron in the state n , and express it 
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by by , then k^r-nk^ , and the criterion becomes k^>Z/a o f which is the same 
as the criterion for excitation of the ground state. 

For ionization of the #th state a-naJZ is obtained. The validity criterion 
is then kJP Z/a o . 

A connection exists between the Bom and Bethe approximations, The Bethe 

approximation (Bethe 1930, Landau and Lifshitz 1958) consists of the assumption 
2 

that 4E , where &E is the energy transfer by the incident electron. With 

this assumption, the cross section can be written 

Q~k* [Ain k* + B], (54) 

with A and B being constants dependent upon the atomic parameters and excitation 
energy for r,. given process. For n*n ! excitation, ^E =(/ / n * -^ / ^,2 )Z 2 /a 2 ^ } and for 


ionization a E = Z Jn a Q +e (cf. 11), Since in ionization dQ/de has appreciable 
value only when e is much smaller or comparable to Z 2 Jn 2 a 2 Q , the assumption 
&E <k 2 is equivalent to > Z/na Q , which is the same as the criterion for 
the validity of the Born approximation in ionization. This means that in the energy 
region where Q , as calculated by the Bom approximation, agrees with Q as given by 
(54), the agreement can be taken as a test for the validity of the Bora approximation. 
In the next section, it will be shown that for ionization of a level with the principal 
quantum number n , the value of k^ at which the Born approximation can be 
approximated by (54) is indeed inversely proportional to n , as expected. 
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4. RESULTS AND DISCUSSION 

To calculate the ionization cross section, the explicit form of T (nlrn, k) 
given by (21) and (28) is substituted into (8) and integration with respect to q and 
the polar angle of k is carried out numerically. No integration with 
respect to the angle of azimuth of if is necessary. With dQ/de determined, a 
final numerical integration is performed according to (10) to determine the 
total cross section Q. 

A quantity of interest is the average energy of an ejected electron for a 
given incident energy. This is given by 



In the following figures and Tables 2 through 6, dQ/de, Q, and e are 

given. 

In Figures 1 through 4 the total ionization cross section for n = 2 through 
5 and all possible / values is plotted as a function of the incident energy. In 
Figure 5, the total cross section is given as a function of energy for n = 1 
through 5, each case being averaged with respect to the azimuthal quantum 
numbers. Figures 6 through 10 are similar to Figures 1 through 5 except that 
C Ejl) Q is plotted as a function of log (E/I) , with EfI the incident energy in 

threshold units, and the abscissa extended to much higher incident energies. If, 




as the incident energy increases, Q , as given by the Born approximation, approaches 
the Q given by (54), then the curves given by Figures 6 through 10 should approach 
straight lines with slopes given by A and with intersections with the ordinates 
given by B . This is noted to be the case. As indicated in Section 3, the agreement 
between the two Q values can be taken as a measure of the validity of the plane wave 
assumption for the incident electron in the Born approximation. Some deviation 
from the straight lines noted in the figures is due to the inaccuracy in the 
numerical integration for the cross section at high incident energies. 

Estimates of A and B for different cases are useful for practical plasma 
and astrophysical calculations. Kingston (1965) gives accurate numerical values 
for A for n = lthrough 5 and all possible / values. Calculation of B is much more 
difficult, and should be estimated from the graphs. 

Studies of the Figures 1 through 10 show that for a given n and incident 
energy, the total cross section tends to increase with increasing / for low and 
intermediate energies, but the trend is reversed at the higher energies. The 
latter situation has been predicted by Bethe (1930). 

These figures also show that similar situations exist with respect to n , 
namely that for low an intermediate energies, for a fixed energy, the cross 
section increases with an increase inw; but at high energies, this quantity 
decreases with an increase in n . The latter situation is consistent with the 
correspondence principle according to which as n increases, the quantal and 
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classical results should agree with each other. Classically, at high incident 
energies, the cross section is inversely proportional to the incident energy, 
and as n increases, A [in (54)] must tend to zero in order to bring Q [as given in 
(54)] into agreement with the classical cross section. The decrease in A as n 
increases is clearly seen in Figure 10. There exist a wealth of experimental 
and theoretical results for the ground state ionization. The subject has been 
recently reviewed by Rudge (1968). In contrast, no measurement is available 
for the excited state ionization, and, with few exceptions, the only available 
calculation other than the Born approximation is the classical calculation. 


Mentioned here is the classical impulse approximation which is valid at 
high incident energies. The original calculation in this method was done by 
Thomson (1912), and recently the method has been revived by Gryzinski (1958) , 
Among many calculations by this method, the one with the least number of 
approximations is discussed here. Gerjoy (1966) has formulated the problem 
of energy transfer between two moving charged particles. Garcia (1969) has 
applied this formulation to the ionization of hydrogenic atoms, where quantum 
mechanical momentum distribution is assumed for the bound electron. He finds 
the differential cross section for a given n, averaged with respect / , to be given 
by 


de 


8e 4 

3 ^/ 3 / 3 



+ 4) (tan' 1 y 


+ 


1 


y 

+ y 2 


) + 


2y(J-4)) 
(1 + y 2 ) 2 } 






1 




/ = 1 + e// , fi t . Ell , y ~(/3 g - ■/), 


(56) 


with e being the electronic charge, and E and e being the incident and the ejected 
electron energies in rydbergs. / = Z 2 e 2 /n 2 a 0 is the ionization potential in By 
of the level n concerned. The total cross section as in (10) is given by 

•e. 


dr = 


f 


max 


dQ 


de 


± de 


(57) 


In Figure 10 (E/I) Q c is plotted as a function of log (Ell) for comparison 

with the quantal cross sections. It is seen that as n increases the quantal calcula- 
tion approaches the classical one. Since both the Born and impulse approximation 
are valid at high energies, this constitutes a verification of the correspondence 
principle. 

As the incident energy E increases, the classical cross section varies as 
E i . This is seen by noticing that 

4 


dQ c 4Tie 
de 


3 ' 2 + 3 J} ’ 


(3,» 1, /«P, 


(58) 


Thus it follows that 


Q, 


i 


dQ c 
Sdj 

dj 3 rp t 


20n e 4 20 n 4 it a 2 
o 

T, ~ 


3 Z 4 E ’ 


(59) 


where only the leading term in the inverse powers of /3 is retained. 
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Equation (59) is in agreement with the asymptotic form of the other classical 
impulse approximation formulas. 

Kingston (1968) has made a classical calculation equivalent to the classical 
calculation contained in (56) and (57). He has made a comparison between his 
results and the privately communicated Born calculation of this author. His 
views are consistent with those given here. However, his calculation is not 
extended beyond = 5. 

In Figure 11 the differential cross section for the ground state ionization 
and incident energy of 25 threshold units is plotted as a function of the energy of 
the ejected electron. The classical expression given by (56) is also plotted for 
comparison. For £// » 1 agreement is noted between the quantal and classical 
calculations. The reason is that formulas for both theories have the same asymptotic 
forms. This can be seen in the following. From (58) it can be seen that for 



» cjl » 1 , 

dQlde ~ (Aire'lEe * ) ( 1 - 2 I/3e ) , E»e» I. (60) 

Omidvar (1969) has similarly shown that the first two leading terms in the ex- 
pansion of dQjde as inverse powers of e/I calculated by the Born approximation 
are identical to (60). 
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Inokuti and Vriens in their classical binary encounter theory find the same 
expression as (60) (private communication, 1968). 

The main disagreement between dQ/de and dQJde is for e// « l where 
the two quantities have different functional forms. Obviously the classical ex- 
pression is not valid in this case. The importance of the region e « / is the fact 
that most electrons in a high energy collision are ejected with small energies. 

The large contribution to the total cross section from the region e « I accounts 
for the difference between the quantal and classical total cross sections. Position 
of the average energy, <?, is shown in Figure 11 by means of a vertical line. 

This position also indicates that large numbers of ejected electrons have small 
energies. 


At this point it is also convenient to discuss the exchange effect. This effect 
is important when e is comparable to E. Contribution of this effect to dQ/de 
is of the order of \/Ee % (Landau and Lifshitz 1958). Integration of l/£e 2 with 
respect to e in the region that e is comparable to E shows that the ratio of 
the contribution of the exchange effect to the total cross section is of the order 
of 1/E , and in the region of the validity of the Born approximation can be neglected. 

In Figures 12 through 14, the differential cross section for the 2s and the 
2p ionization for the incident energies of 4, 25, and 100 threshold units are 
plotted versus energy of the ejected electrons. The positions of the average 
energy of the ejected electrons, e", [cf. (55)], are indicated by vertical lines. 
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Concerning the accuracy of the numerical results presented here* two 
checks have been made to rule out a systematic error in the calculations. One 
is the comparison between the present and previous results (Omidvar 1905) for 
a given n and n * In the present calculation, the % -axis is along k, while in 
the old calculation this axis is along the momentum transfer. When the two 
results have been summed with respect to m, agreement between them is found. 
The next check has shown agreement between two sets of results where, in one 
set, summation has been performed with respect to « , and in the second set 
with respect to / . 


Gaussian quadratures have been used as a method of integration, and for 
each Q as a function of the parameters kp n, and /; the number of intervals 
of the integrations has been increased until successive values of the Q have 
fallen within 1% of each other. An exception to this rule is made for n~ 5 
where the differences are 2 to 3%. 

2 

Since the range of integration increases as kj increases, numerical 

2 2 
values at low k ^ are more accurate than those for large kj , . 


5. SPECTRUM OF THE SECONDARY ELECTRONS IN THE PASSAGE OF AN 
ENERGETIC ELECTRON IN AN ATOMIC HYDROGEN GAS 


In this section are considered the number of ejected electrons per unit of 
their energy range, and the total number of ejected electrons, when an energetic 
electron enters an atomic hydrogen gas and is brought to rest through ionization. 
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We assume that the atoms are in their ground state. The results obtained can 
be applied to the case of hydrogen gas if the few electron volts binding energy of 
the hydrogen molecules can be neglected compared to the energy of the primary. 

If is defined as the energy lost by the primary during the course of a 
single collision due to all inelastic processes, in anology to (55) it can be 
calculated by the expression 


00 

r Q max 
/ (/ + <0 

I. 

jQin * 

(61) 

n = 2 

*0 

J 

! 



where Q f} is the cross section for excitation of the n^ 1 level, E 0 and E n are 
energies of the ground and n*- 1 levels, / is the ionization potential, and Q( n 
the total inelastic cross section. If (IN is the number of ions produced when a 
primary of energy E loses an amount of energy clE , we must have dN - (dEf&E) ( QlQjn ) 
with (2, as before, being the total ionization cross section. The number of 
ejected electrons having energies between e and e + de is a fraction 
( dQ/de ) dejQ of dN. This number can then be written 


■# 


* 


dN dtdE- * 2 !* -J- de dE 


dedE 


Q Q in *E 


dQ/de 
Qin AE 


dedE . 


(62) 


If the initial energy of the primary is . the number of electrons with 

energies between € and e + de is given by 
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1 + e 


cte Q ln *E 


For evaluation of (63), a knowledge of Q^E is necessary. Accurate 
evaluation of this quantity is difficult. However, at high incident energy and 
within the approximation applied to the derivation of (54), it can be written 
(bethe 1930, Landau and Lifshitz 1958) 


Q ln A E ~ m/E) in C2E/J) , 


with Q. n in units of ira* and E in rydberg. In (64) Z is the number of atomic 
electrons, and 7 is the average excitation energy in rydberg. For the ground 
state atomic hydrogen 7-1.1 05 Ry . 


Although the Born approximation is inaccurate at low incident energies , it 
is likely that the ratio of two quantities calculated according to the Born approxi- 
mation will be more accurate than the quantities by themselves. The ratio in 
question is Q/Q.^E given by (62). Making use of (54) and (64), we can write 


b e “ 8 Z 


dQjde Q4 to £•+•/?)</£ 
Q x In (IE ID 


/ + € 


with / the ionization potential of the bound electron. Using (65), dN/de for the 
ground state atomic hydrogen has been evaluated numerically. Figure 15 is a 
plot of dN/de versus e for various initial primary energy E 


The total number of secondary electrons for a given E is obtained by 

ftTCtX 

integrating (65) with respect to e . This is given by 
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A In E + B 

1K2W dE 




w|. £( (L (1L 


where E,(x) is the exponential integral. 


Using (66), a plot of N for the ground state of atomic hydrogen is given in 
Figure 16. Dalgarno and Griffing (1958) have calculated N due to ionization by 
the primary and all higher generation electrons. This result is also shown in Figure 16. 
As the incident energy increases, the value of N given by Dalgarno and Griffing in- 
creases more rapidly than the value of N given by (66), mainly due to the ionization 
produced by the higher generation electrons . 


The energy spectrum given here can also be used to deduce the energy 
spectrum of the tertiary electrons due to energetic secondary electrons, etc. 


It is of interest to know that when e is several threshold units, for any target 
atom or molecule dQ/de to a good approximation is given by the Rutherford formula 
which is given by the leading term of (60) multiplied by Z . Making use of (54) and 
the leading term of (60) in (65) results in the following analytic formula: 


bN 1 f E m 

be 2e 2 / 

*7+6 


ln(2E/D 


=4-7 k f 2 J f*) 7 $ f > / 
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Table 1 


Atomic dipole moments in units of eY for Is —’tip transistions. 
23456789 10 


P 0.7449 0.2983 0.1759 0.1205 0.0896 0.0701 0.0569 0.0474 0.0403 



l. 



Table 2. 


Ionization of the n»2 levels, kj, is the ratio of the speed of the incident electron 
to the speed of the bound electron, k 2 gives the incident energy in threshold units. 
Q is the total cross section in units of na* , and e is the average energy of the 
ejected electron in Ry . / is the azimuthal quantum number. For each / the 
relevant quantities have been averaged with respect to the magnetic quantum 
number m , 


U 0 /« / 


,v r 

n 4 Q 

2 

n e 

n 4 Q 

2 

n e 

1.2 

0.668 

0.163 

0.643 

0.155 

1.4 

1.017 

0.321 

1.083 

0.303 

1.6 

1.103 

0.465 

1.234 

0.443 

1.8 

1.074 

0.596 

1.227 

0.572 

2.0 

1.001 

0.714 

1.152 

0.689 

2.5 

0. 792 

0.965 

0.902 

0.936 

3.0 

0.620 

1.165 

0.692 

1.130 

3.5 

0.493 

1.328 

0. 541 

1.285 

4.0 

0.399 

1.462 

0.431 

1.414 

4.5 

0.329 

1.572 

0.352 

1. 523 

5.0 

0.276 

1.663 

0.292 

1.619 

6.0 

0.202 

1.814 

0.210 

1.783 

7.0 

0. 155 

1.911 

0.159 

1.923 

8.0 

0.122 

1.969 

0.124 

2.043 

9.0 

0.099 

1.999 

0.100 

2.145 

10.0 

0.082 

2.009 

0.082 

2.229 

15,0 

0.040 

2.142 

0.039 

2,496 

20.0 

0.023 

2.120 

0.022 

2.575 

25.0 

0.015 

2.117 

0.014 

2. 59 

30.0 

0.011 

2.125 

0.010 

2.59 
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Table 3 


Ionization of then » 3 levels. 

Notations and explanations are the same as in Table 2. 


f « 0 / = 1 l = 2 


k T 

-4 

n Q 

V 

2. 

n e 

-4 _ 

n Q 

2- 

n e 

-4 

n Q 

2- 

n e 

1.2 

0.708 

0.167 

0.701 

0.161 

0.665 

0.150 

1.4 

1.034 

0.327 

1.072 

0.317 

1.149 

0.297 

1.6 

1.101 

0.473 

1.161 

0.462 

1.307 

0.437 

1.8 

1.059 

0.607 

1.127 

0.594 

1.290 

0.560 

2.0 

0.978 

0.729 

1.046 

0. 714 

1.199 

0.689 

2.5 

0.760 

0.990 

0.816 

0.971 

0.917 

0.944 

3.0 

0.587 

1.200 

0.630 

1.178 

0.692 

1.146 

3.5 

0.463 

1.372 

0.494 

1.349 

0.533 

1.311 

4.0 

0.373 

1. 514 

0.396 

1.495 

0.421 

1.449 

4.5 

0.306 

1.628 

0.324 

1.625 

0.340 

1. 569 

5.0 

0.255 

1.717 

0.270 

1.743 

0.281 

1.674 

6.0 

0.185 

1.752 

0. 195 

1.937 

0. 200 

1.853 

7.0 

0. 140 

1.792 

0. 148 

2. 064 

0.150 

2. 019 

8.0 

0.110 

1.833 

0.116 

2.162 

0.116 

2.145 

9.0 

0.089 

1.883 

0. 093 

2.202 

0.093 

2.270 

10.0 

0.074 

1.974 

0.077 

2.239 

0.076 

2.389 

15.0 

0.035 

2.218 

0.036 

2.361 

0. 035 

2.771 

20.0 

0.021 

2.365 

0. 021 

2.510 

0.020 

2. 963 

25.0 

0.014 

2.382 

0.014 

2.609 

0. 013 

3.045 

30.0 

0.010 

2.388 

0.010 

2.614 

0. 009 

3.139 
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Table 4 


Ionization of the n - 4 levels. 

Notations and explanations are the same as in Table 2. 


T -4- 2- -4^ 

n Q n e n Q 

1.2 0.726 0.168 0.724 

1.4 1.048 0.328 1.070 

1.6 1.108 0.476 1.143 

1.8 1.059 0.611 1.099 

2.0 0.973 0.735 1.013 

2.5 0.749 1.001 0.783 

3.0 0.575 1.216 0.602 

3.5 0.450 1.387 0.472 

4.0 0.360 1.515 0.378 

4.5 0.294 1.602 0.309 

5.0 0.244 1.660 0.258 

6.0 0.175 1.707 0.186 

7.0 0.133 1.844 0.141 

8.0 0.105 1.954 0.110 

9.0 0.085 2.056 0.089 

10.0 0.070 2.118 0.073 

15.0 0.034 2.422 0.034 


1 

l 

as 2 

l 

= 3 

2 - 

n € 

-4 

n Q 

a*« 

"4 ~ 

n Q 

2 - 

n e 

0.164 

0.709 

0.157 

0.661 

0.144 

0.323 

1.116 

0.311 

1.186 

0.290 

0.469 

1.219 

0.454 

1.358 

0.431 

0.603 

1.185 

0.585 

1.335 

0.564 

0.726 

1, 098 

0.706 

1.233 

0.686 

0.989 

0.849 

0.965 

0.930 

0.946 

1.203 

0.649 

1.174 

0.694 

1.154 

1.383 

0.504 

1.344 

0. 530 

1.325 

1.543 

0.401 

1.484 

0.416 

1.470 

1.690 

0. 326 

1.608 

0.335 

1.595 

1.820 

0.270 

1.725 

0.274 

1.705 

1.966 

0. 194 

1.956 

0.194 

1.898 

2. 015 

0.146 

2.166 

0.144 

2. 078 

2.044 

0.113 

2.274 

0.112 

2.214 

2.099 

0.091 

2.351 

0.089 

2.360 

2.189 

0.074 

2.417 

0.072 

2.496 

2.481 

0.035 

2.615 

0.033 

2.977 
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Ionization of the n = 5 levels. 

Notations and explanations are the same as in Table 2. 
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FIGURE CAPTIONS 


Figure 1. Total cross section for the ionization of the n * 2 level of atomic hydrogen, 
Cross sections, in units of n A v(? • are averaged with respect to the magnetic quantum 
number m for each azimuthal quantum number / * For low and intermediate incident 
energies Q for / * X is larger than Q for l = 0 for a fixed energy, but at high energies 
it falls below / - 0 values (cf, Figure 6), 

Figure 2, Total cross section for the ionization of the n = 3 level. Similar to the 
n = 2 case, for low and intermediate energies Q increases as / increases, The 
situation is reversed at high energies , 

Figure 3, Total cross section for the ionization of the n » 4 level. 

Figure 4. Total cross section for the ionization of the n - 5 level. 

% 

Figure 5. Total cross section for the ionization of the n ” 1, 2, 3, 4, and 5. For 
each/? , the cross section has been averaged with respect to the azimuthal quantum 
number / , E is the energy of the incident electron, and I is the ionization 
potential for the level concerned. With the scaling shown in the Figure, and for 
a given low energy, the cross section increases with the increasing order of n , 

The situation reverses for snergies higher than a few threshold units. 

Figure 6. A plot of (E/I) Q versus log l0 (E/f) for ionization of the n - 2 

and l = 0, 1 states. At high energy the graphs approach straight lines. Slopes 


and intersections of these lines with the ordinate give the two parameters 
necessary for the high-energy calculation of Q. 

Figure 7. A plot of (/i/7) Q versus log iQ (E/I) for ionization of the n « 3, 
l e 0, 1, 2 states. Notations are the same as in Figure 6, 

Figure 8. A plot of (EfJ) Q versus log iQ (E/I) for ionization of the n s 4 states. 
Notations are the same as in Figure 6. 

Figure 9. A plot of (E/I) Q versus log i0 (E/I) for ionization of the n ** 5 states. 
Notations are the same as in Figure 6. 

Figure 10. A plot of (E/I) Q versus log lQ (E/I) for n « 1, 2, 3, 4, and 5 
levels. For each m Q has been averaged with respect to the different / . The 
classical formulais that of the impulse approximation with quantum mechanical 
momentum distribution assumed for the bound electron. As n increases the 
quantum mechanical results approach those given classically in accordance with 
the correspondence principle. 

Figure 11. A comparison of the results of the Born, the asymptatic Born, and 
the classical theories for the differential cross section with ground state ionization 
of hydrogen. The asymptatic Born is valid for e > / . The incident energy is 
25 rydbergs. The average ejected electron energy ? is shown by a vertical line. 


Figure 12. Differential cross section for the n - 2 level ionization for an 
incident energy of 4 threshold units. Positions of 7 (2s) and e(2p) ax’e shown 
by vertical lines . 

Figure 13. Differential cross section for n ~ 2 and an incident energy of 25 
threshold units. Notations are the same as in Figure 12. 

Figure 14. Differential cross section for n » 2 and an incident energy of 100 
threshold units. Notations are the same as in Figure 12. 

Figure 15. Number of the secondary electrons per unit range of their energy, 
dN/ 3e, versus e for different primary energies. The primary energies 
through E are 10, 20, 30, 40, 60, 80, and 100 units of the ionization potential, 
respectively. 

Figure 16. Total number of the ejected electrons versus the primary energy. 
Curve 1 is the result of the present calculation. Curve 2 is due to Dalgarno 
and Griffing (1958). 
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